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We study the motion of neutral test particles in the gravitational field of two charged black holes 
described by the extreme Reissner-Nordstr0m dihole metric where the masses and charges of the 
black holes are chosen such that the gravitational attraction is compensated by the electrostatic 
repulsion. We investigate circular orbits in the equatorial plane between the two black holes with 
equal masses as well as the case of circular orbits outside this symmetry plane. We show that the 
first case reduces to an effective two-body problem with a behavior similar to a system described 
by the Reissner-Nordstr0m spacetime. The main focus is directed to the second case with circular 
orbits outside the equatorial plane. 
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I. INTRODUCTION 

The study of binary black hole systems is of great im- 
portance in the context of the emission of gravitational 
waves. The direct proof of the existence of such waves is 
still lacking and is one of the central goals in modern as- 
trophysics. Their discovery would result in another great 
piece of support for the validity of Einstein's general rel- 
ativity. Indirect evidence of gravitational waves has al- 
ready been obtained in radio observations of binary pul- 
sars, such as the famous binary pulsar PSR 1913+16, for 
which Hulse and Taylor won the Nobel Prize in physics; 
see PQ. Furthermore, pairs of black holes could play an 
important part in the evolution of galaxies and the whole 
Universe; see for example Pallerola [2]- 

This motivates a methodical analysis of multi black 
hole systems. However, for such systems no general ana- 
lytic expression exists, in contrast to a single black hole. 
In addition, because of their character, these systems 
have to be described by nonstatic spacetimes. The ef- 
forts of many authors to find numerical solutions of Ein- 
stein's field equations also shows the relevance of this 
issue. There are, for example, approximated metrics for 
binary black hole systems in a circular orbit with widely 
separated, nonrotating black holes by Alvi [3] or a global 
approach by Gourgoulhon et al. [4]. 

The extreme Reissner-Nordstr0m (RN) dihole metric 
describes a system of two charged black holes with com- 
pensating gravitational and electrostatic forces and is a 
version of the two-center problem in general relativity. 
Thus, a static and analytic spacetime for this special case 
of a double black hole system exists, which can be used as 
a comparatively simple model metric to analyze multi- 
black hole systems. 

A first step in the study of curved spacetimes is to con- 
sider the motion of test particles. The large diversity of 
null geodesies in the meridian plane was already pointed 
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out by Chandrasekhar [5] . In contrast to a corresponding 
Newtonian system, the relativistic motion can be chaotic 
as was shown by Yurtsever [BJ. Here, we are particularly 
interested in the existence of lightlike and timelike circu- 
lar geodesies. Because of their simplicity, they provide 
a good starting point for exploring a metric. In general, 
they play an important role in analyzing a spacetime, 
as indicated by numerous publications on this topic, in- 
cluding Pugliese et al. [7j, Chowdhury et al. [S], or Bini 
et al. 0. 

The compensating attractive and repelling forces in 
the dihole system are motivated by classical mechanics. 
However, a strict general relativistic calculation shows 
the consistency of this picture with the Einstein-Maxwell 
equations; see 0[1O]. Thus, in this context, it is also of 
interest how far the dynamics of the classical and the 
general relativistic systems differ. In this paper, we ana- 
lyze circular orbits in the general relativistic case. A brief 
overview of the classical analogon for massive particles is 
given in Appendix [B] and shows a remarkable relation- 
ship to its general relativistic counterpart, but also some 
significant differences. 



For discussing the dynamics in the extreme RN di- 
hole spacetime, we distinguish between two cases. The 
first one is the motion in the equatorial plane of the ex- 
treme RN dihole metric for equal masses. In this case, 
we have an effective two-body problem similar to the case 
of a naked singularity of the simple Reissner-Nordstr0m 
spacetime, which we will show with calculations related 
to the paper by Pugliese et al. [TT]. They studied the 
existence of circular orbits in the RN metric as well as 
their stability for black holes and naked singularities by 
discussing the effective potential. They also investigated 
the motion of a charged particle in the RN metric |12) . 
Our main contribution is the discussion and analysis of 
the second case, namely circular orbits outside the equa- 
torial plane, which are, to the best of our knowledge, not 
covered in previous literature. 
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II. EXTREME REISSNER-NORDSTR0M 
DIHOLE METRIC 

The extreme RN dihole metric reads 



ds 2 = g^dx^dx" — 



df 
U 2 



+ U 2 {dx 2 + dy 2 +dz 2 ) (1) 



with 



U(x,y,z) = 1 + 



Mi 

-y/a; 2 + y 1 + (z - l) 2 

M 2 

^/z 2 + y 2 + (z + l) 2 



(2) 



(see [5]) and describes a system of two extreme 
RN black holes with masses Mi and M2 at positions 
ri = (0, 0, z\ — 1) and r 2 = (0, 0, z 2 = —1) and charges of 
the same sign. The gravitational attraction is compen- 
sated by the electrostatic repulsion. For the form and 
properties of the RN metric, see for example [TUl EH [H] . 
The metric is given in geometric units, where the 
speed of light c and Newton's gravitational constant G 
are normalized to unity. The extreme RN dihole metric 
is a special case of the Majumdar-Papapetrou spacetimes 
(see for example (10) V which describe an arbitrary num- 
ber of extreme RN black holes with compensating gravi- 
tational and electrostatic forces. This interpretation goes 
back to Hartle and Hawking [T5] . 

Since the system is axisymmetric with respect to 
the z axis, we use cylindrical coordinates. The metric 
then transforms to 

ds 2 = ~ + U 2 (dp 2 + p 2 dip 2 + dz 2 ) (3) 
U 4 

with 

vV + (z - l) 2 vV + tz + l) 2 

We call z — the equatorial plane. A plane containing 
both singularities is denoted a meridian plane. 

From Eq. ^j, we obtain the Lagrangian of a neutral 
test particle 



U' 



+ U 2 (p 2 +p 2 <p 2 + z 2 ) 



(5) 



Here, the dot denotes differentiation with respect to an 
affine parameter A, which in the case of a timelike particle 
can be identified with its proper time. Additionally, the 
constraint g llu x tl x 1 ' = k has to be fulfilled, with k = 
for null geodesies and k = — 1 for timelike geodesies. In 
this paper, the metric signature sign(g M „) = +2 is used. 

Combining this constraint and Eq. ^ immediately 
yields £ = k/2 = const. Further constants of motion 
follow from the cyclic coordinates t and tp: 



E 



d£ 
~dt 



t 

U 2 ' 



dip 



p 2 U 2 <p. 



(6) 



E represents the energy of a particle and L z its angu- 
lar momentum with respect to the z axis. With these 
constants, Eq. ^ can be rewritten as 

1 p 2 + h 2 + V cfi ( Pl z) = 1 ^- (7) 



T 2 
with the effective potential 

V cS (p,z) 



2\pzm(p,z) u 2 { Pl z)! (S| 

and U(p,z) from Eq. Q. We choose the initial direc- 
tion of a particle with respect to a local tetrad. For the 
metric ([3]), the tetrad vectors can be chosen as 

1 

U' 
1 

U' 



e (t) - Ud t , e (p) - -,.„. 



1 



-A 



(9) 



Then, an initial direction u can be given by 

u = ±w (t) e (4 ) + tpn (10) 

with n — sin x cos £ e( p ) +sinx sin£ e^) +cos \ e (z)- The 
sign in front of vA determines the direction of time. For 
a lightlike particle, we have = ip = 1, and for a time- 
like particle u^> = 7 and ip — /?7, with the local velocity 
/3 and the corresponding Lorentz factor 7 — (1 — fi 2 )^ 1 / 2 . 
In the case of a timelike particle, u can be identified with 
its four-velocity. 



With respect to the initial direction ( 10 ) the constants 
of motion E and L 7 can also be written as 



E = 



± 
U' 



L z = ip sin x sin £ pU. 



(11) 



It is a priori obvious that there is a point of equilib- 
rium z equ on the z axis, where the gravitational forces 
of the two masses Mi and M 2 on a timelike test parti- 
cle (k = —1) compensate each other. A timelike particle 
with vanishing initial velocity at this point remains at 
rest, which is impossible for a lightlike particle. We can 
calculate z cqu from Eq. ([7]) by considering a test particle 
with initial conditions p{6) = p(0) = 0. In this case, L z 
from Eq. (11) vanishes and thus we have 



ii 2 + V eS (z) = ^- with 1^(3) = _L_ (12) 



and 



U{z) = 1 + 



Mi M 2 



1-z 1+z 



(13) 



For null geodesies, the effective potential would vanish. 
Figure [l] shows V e g(z) for some combinations of the two 
masses Mi and M 2 . The point of equilibrium z cqu is 
given by the local maximum in the interval (— 1;+1). 
The condition dV c ff/dz — yields 



^cqu 



(<?) = 




5^1 
9 = 1 



(14) 
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FIG. 1. (color online) Effective potential V e ff(z) from Eq. ( 12 1 
for three different combinations of the mass parameters Mi 
and M2. The local maximum is always in the interval 
(— 1; +1) and marks the point of equilibrium z cqu . 




FIG. 2. (color online) Point of equilibrium z c qu ( q) from 
Eq. (141 and the singularity z s i ng (q) from Eq. (34 1. It is 

z cq u(q = 0) = +1 and z equ (q ->• 00) = -1 



mass is than the other, the closer z c 



The bigger one 
lies to that mass. The 



quantity z s ing has exactly opposite properties and becomes 
important in Sec. |IV| 



with the mass ratio q := M1/M2. The value z oqu depends 
only on the ratio of the two masses, and is shown in Fig. [2] 



III. CIRCULAR ORBITS IN THE EQUATORIAL 
PLANE 



Since the extreme RN dihole spacetime is axisymmet- 
ric, we define a circular orbit as an axisymmetric closed 
orbit. Therefore, all circular geodesies fulfill dV c s/dz = 
and dV cS /dp = with V eS (p,z) from Eq. The first 
condition leads to dll/dz — 0, both for null and timelike 
geodesies. This results in the expression 



Mi 
W 2 



l + z 



1 - 



42 



p 2 + (z + l) 2 



3/2 



(15) 



which shows that circular orbits can only exist in the 
range z e (—1; +1), because both masses have to be pos- 
itive. Furthermore, for z = 0, Eq. ( [15] ) becomes indepen- 
dent of the radius p, and it follows Mi = M2. Therefore, 
a particle in the equatorial plane (z — 0) with zero initial 
velocity components normal to this plane remains there, 
if the masses Mi and M 2 are equal. In this case, we 
define M ;=M 1 = M 2 and Eq. Q simplifies to 



1 



P 2 + V cfi (p) = 



2 



with 



and 



VM = 



2 \p 2 U\p) 



K 

uKp) 



U(p) = 1 + 



2M 



(16) 



(17) 



(18) 



A. Null geodesies 



The effective potential (17) for null geodesies is shown 
in Fig. [3] for different values of the mass M and the an- 
gular momentum L z . The radii of circular orbits follow 
from the extremal points of that effective potential. The 
condition dV c ff/dp = for k = leads to the expression 

(p 2 + l) 3 / 2 =2M(p 2 -l), (19) 

which has solutions only for p > 1. Substituting 
i] := p 2 + 



1 in Eq. ( 19 ) yields the cubic equation 
rf - AM 2 if + 16M 2 r? - 16M 2 = 



(20) 



with the constraint rj > 2 because of p > 1. This cubic 
equation can be solved using Cardano's formulas. Taking 
77 > 2 into account, we obtain 



Vi(M) = - 



with 



M 2 + 2M^M 2 - 3 cos 
M 2 - 2M^M 2 - 3 cos 



$(M) := arccos 



27 - 36M 2 + 8M 4 



8M V / (M 2 - 3) 3 



(21a) 
(21b) 

(22) 



The quantities r\\ and 772 are defined only for M > M 
with 

M ■= \ — ~ 1.83712, (23) 
V 8 

which follows from the roots of the discriminant of the 



cubic ( 20 1 . The two corresponding photon orbits with 
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FIG. 3. (color online) Effective potentials for a lightlike particle with different angular momenta L z for three different values 
of M. (a) M = 1: There are no extremal points, (b) M = M » ^27/8: There is one photon orbit, (c) M = 2.5: There are 
two photon orbits. The radius p„ with V e s{p*) = V e ff(p p i) and p» 7^ p p i is also shown. 



radii p p \ and p P 2 then follow from p = y/rj — 1. The 
radius p p i(M) is a local maximum of T4g(p) and thus the 
orbit is an unstable circular orbit. The radius p P 2(M) is 
a local minimum of V c ff(p) and thus the orbit is stable in 
a local sense along the p direction. Both photon orbits 
are shown in Fig. |4j 
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FIG. 4. (color online) Photon orbits p p i/2(M) in the equa- 
torial plane. They do not exist for M < M and coincide for 
M = M. 



The existence of a local maximum in the effective po- 
tential implies the existence of an asymptotic trajectory 
to this photon orbit. To see this we consider an observer 
at distance po in the equatorial plane that emits a light 
ray at angle £ (see Fig. [5]). Such an asymptotic orbit has 
to fulfill 



Kiax := Ves(ppl) 



2 



(24) 



[see also Fig.[3jc)] . The constants of motion E and L z can 
be expressed with respect to the local tetrad according 
to Eq. ( 11 1 with z = and \ = tt/2. Equation (24) then 
reads 



p 2 a U 2 (Po) sin 2 £ 
2^i^ 4 (P P i) 
and from this, we obtain 



2U*(p ) 



£(po) 



arcsin ( 



7r — arcsin 



PpiC/ 2 (p P i) 
PoU 2 (p ) 



Pa > P P i 



( Ppl u 2 ( Ppl ) \ . < 



(25) 



(26) 



with U(p) from Eq. (18). Typical curves of £(po) are 
shown in Fig. [6j The function is only defined if the pho- 
ton orbit p p i exists. For an observer at this orbit, the 
angle is £ = 90° in accordance with the similar situation 
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FIG. 5. (color online) A light ray is emitted in the equatorial 
plane from po at an angle £ and asymptotically approaches 
the photon orbit p p \. 



FIG. 8. (color online) Effective potentials for a timelike par- 
ticle in the case M — 1.5 with different angular momenta 
L z . In this and the following figures, the radii pi sco and pi uco 
stand for the last radially stable and unstable circular orbits, 
and are obtained from Eq. ( 29 1 . Further information on their 
nomenclature is in the text. 
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FIG. 6. (color online) The angle £(po) from Eq. ( 26 1 in degrees 
for M — M — yW/8 and M = 2.5. 



FIG. 9. (color online) Effective potentials for a timelike par- 
ticle in the case M = M = y/27/8 « 1.837 with different 
angular momenta L z . 
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FIG. 7. (color online) Effective potentials for a timelike par- 
ticle in the case M = 1 with different angular momenta L z . 
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FIG. 10. (color online) Effective potentials for a timelike par- 
ticle in the case M — 2.5 with different angular momenta 
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FIG. 11. (color online) The functions E 2 (p)/2 and L 2 z {p)/2 
for M = 1 as an example of the interval < M < M* « 1.029 
see Eq. p0| ]. For all L^/2, there is one circular orbit (see 
also Fig.[7|. For M = M*, there would be a reflection point in 
the diagrams. In this and the following figures, the stability 
refers to the radial direction. 



FIG. 12. (color online) The functions E 2 (p)/2 and L 2 z (p)[2 
for M = 1.5 as an example of the interval M, < M < M. 
There can be one, two, or three circular orbits. The area for 
L 2 /2 with three circular orbits is bounded (see also Fig. 
These boundaries are given by the radii pi ECO and pi uco , which 
follow from solving Eq. ( 29 1 . 



in the Schwarzschild metric (compare with Miiller |16j). 
which can also be easily expanded to the RN metric. 



B. Timelike geodesies 

The circular orbits of a timelike test particle follow 
from the extremal points of the effective potential in 
Eq. ( IT I for k = — 1. This function is shown in Figs. [7 10 



for different masses M and angular momenta L z . 

Here, the position of circular orbits depends on L z in 
contrast to the lightlike case. In the potential, L z ap- 
pears only quadratic; thus, when considering L 2 , we do 
not have to distinguish between the circumferential direc- 
tion on the circular geodesic. For the search of extremal 
points, we have to solve the equation dV c g/dp = 0, which 
can be resolved with respect to L 2 /2. Again, we substi- 



tute 77 



1 and obtain 



1,2(7?) M [if + AMif/ 2 + (AM 2 - 2)tj 2 - 8r] 3 / 2 ] 



77(773/2 - 2Af?7 + AM) 
M [(1 - 8M 2 )?7 + AMrj 1 / 2 + AM 2 } 
7;(7? 3 /2 - 2Mrj + AM) 



(27) 



Next, we replace L 2 /2 in the effective potential from 
Eq. (17). Using p = 0, we obtain from Eq. (16) the 



expression 



E 2 ( V ) 



r]{rf /2 + 2M) 



2(7^/2 + 2M) 2 (rf/ 2 - 2Mt] + AM) ' 



(28) 



These two functions give the angular momentum and 
the energy of a particle on a circular or bit wit h radius 
p = yjr\ — 1. They are plotted in Figs. |Hr[l4| as func- 



They are plotted in Figs. 11 
tions of p for different masses M. The functions have 
singularities at the photon orbits p p i/2 and are negative 
in between. The angular momentum has the root 77 = 1, 
which corresponds to the radius p — 0. This value rep- 
resents the circular orbit at the center point of the two 
masses with vanishing radius and coincides with the point 



of equilibrium from Eq. ( 14 ) 



From Figs. [lT}{14| it is obvious that there is a different 
number of circular orbits for different values of M and 
L 2 /2. For an exact analysis, the extremal points have to 
be calculated. Both dL 2 /drj = and dE 2 /dr] = lead 
to the condition 

f(rj) := rf -6Mrj 5/2 +3if + 22Mrf /2 + 16M 2 = 0. (29) 

Additionally, as in the lightlike case, we will analyze the 
stability of the circular orbits in the p direction. To do 
so, we consider d 2 V e g/dp 2 = and again use the substi- 
tution 77 = p 2 + 1 and Eq. ([27]). The result is also the 
expression from Eq. (29). Thus, the extremal points of 



L 2 /2 and E /2 are linked directly to the reflection points 
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FIG. 13. (color online) The functions E 2 (p)/2 and L 2 z {p)/2 
for M = M = -s/27/8 w 1.837. There can be one, two, or 
three circular orbits. The area for L 2 /2 with three circular 
orbits no longer has an upper boundary (see also Fig.[9|. The 
singularity is given by the two degenerate photon orbits p p i/2 
and the local minimum by pi sco , which follows from solving 
Eq. ((291. 



FIG. 14. (color online) The functions E 2 (p) /2 and L 2 (p) /2 for 
M — 2.5 as an example of the interval M > M. The behavior 
is qualitatively similar to the case M = M. The correspond- 
ing effective potential is shown in Fig. |10| Between the two 
photon orbits p p i and p P 2, there cannot be any circular orbits. 
The local minimum is given by pi sco and follows from solving 
Eq. p9|. 



of Ves(p), and we obtain the radial stability of a circu- 
lar orbit from the sign of their deri vative (such as those 
shown in the plots of E 2 /2 in Figs. 



In Fig. 15 we present f(rf) of Eq. ( 29 1 for several values 
of M. The existence of roots depends on the explicit 
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14) 



value of M. As shown in Fig. 15 there is a characteristic 
value that divides the M range in an interval with 
no roots of f(rf) and in an interval with two roots. For 
M = M t , there is exactly one root. This mass M* follows 
from f(rj) = and f'(rj) = and is given by 



(13 + \/m)v / 710" 



1.02949. (30) 



Equation ([29]) is discussed numerically. If roots exist 
(M > M*), they are called p\ sco and pi UC o, and are shown 
in Fig. [16] The subscript "lsco" (last stable circular or- 
bit) is motivated by the analogy to the corresponding ra- 
dius known from the Schwarzschild jT7] or the RN space- 
time |llj . Here, however, we refer only to the radial sta- 
bility. Furthermore, pi SC o is the boundary where the ra- 



dial stability character changes (compare with Figs. 12 
14). For p < pisco, the circular orbits are unstable, and 
for p > pisco, they are radially stable. The radius pi uco 
exhibits the opposite behavior, and for this reason we use 
the subscript "luco" (last unstable circular orbit). 

In fact, for p < pi uco , there again appears a domain 



with radially stable orbits, and the innermost stable cir- 
cular orbit is given by rj = 1, which corresponds to p = 0. 
Nevertheless, we will call the local minimum of E 2 /2 and 
Ll/2 the last stable circular orbit. For M > M, the nu- 
merical solution of Eq. ( 29 ) corresponding to pi UC o is be- 
tween the two photon orbits p p i and p P 2; thus pi UC o does 
not exist in this M range. The radii pi SC o and pi uco are 
also marked in Figs. and Figs. [1§{T4| 

The regions with different numbers of circular orbits 
are bounded by pi sco and pi uco or 77^0 and ?7i uco , respec- 
tively Inserting these numerical values into Eq. ( 27 ) , we 
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FIG. 15. (color online) The function f(rj) from Eq. (29 1 for 



different mass values. The number of roots depends on M 
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FIG. 16. (color online) Numerical solution of Eq. (291 
tions exist for M > M« 



Solu- 



For M > M, the dash-dotted branch 
would lie between the two photon orbits p p i/2 (see also Fig. [4]) 
but would have no physical meaning. 




FIG. 18. (color online) Local velocity /3(p) on a circular orbit 
with radius p. This velocity vanishes for p = 0, and becomes 
P = 1 at the photon orbits p v \n- 



After resubstituting, we obtain 



obtain corresponding boundary lines in the (M,L 2 /2) 
plane. These are displayed in Fig. 17 The figure shows 
the number of circular timelike orbits with dependence 
on the parameters M and L 2 /2 of the effective potential 
from Eq. (TP? I. 



To remain on a circular orbit a particular velocity is 
necessary. From Eq. (Ill, we obtain the angular momen- 
tum L z of a timelike particle. We restrict ourselves to cir- 
cular orbits (£ = 7r/2) in the equatorial plane (x = tt/2) 
and obtain 



with T) := p + 1. Comparing this relation with Eq. (27) 



leads to the necessary local velocity fi{rj) on the circular 
orbit with radius p{rj) depending on the parameter M. 
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FIG. 17. (color online) Number of circular timelike orbits with 
dependence on the parameters M and L 2 z /2. The different 
regions are colored uniformly. In parentheses, we indicate 
whether the orbits are radially stable (s) or unstable (u). The 
order is given by their appearance with increasing p in the 
effective potential from Eq. dlTI. 



2Mp 2 



2M + (p 2 + l) 3 / 2 ' 



(32) 



This function is shown for several masses M in Fig. [18] 
At p — 0, we have /3 = 0, which again is the degenerated 
circular orbit at the point of equilibrium. At the photon 
orbits Ppi/2, the velocity is (3 — 1, and in between it 
becomes /3 > 1. Thus, there are no timelike circular 
orbits in the range p P 2 < p < p p i. This agrees with L 2 /2 
and i? 2 /2 from Eqs. (|27"|) and (28), which are singular at 



the photon orbits and undefined in between. 



IV. 



CIRCULAR ORBITS OUTSIDE THE 
EQUATORIAL PLANE 



In Sec. Ill 



(151, 



we have shown, by means of Eq 
that lightlike and timelike circular orbits only exist for 
z 6 (—1; +1). For z = 0, this equation was independent 
of p and implied Mi = M 2 . For < \z\ < 1, Eq. (15) can 
be rearranged to 



P 2 = 



1 



4z 

2/3 



1+2 

1-2 



2/3 



(33) 



which now shows a direct dependence on the radius p 
and the height z for the motion with z = 0. The exis- 
tence of circular orbits outside the equatorial plane can 
be understood classically, at least for massive particles, 
by analyzing the gravitational forces; see Appendix [B| 
Photon orbits are a purely relativistic effect. 

The radius p(z) depends only on the mass ratio 



q := M1/M2, and is shown as a function of z in Fig. 19 
for several values of q. For equal masses (q = 1), all 
heights z with < \z\ < 1 are possible candidates for 
circular orbits. In the case of unequal masses (q =/= 1), 
there is a domain without allowed radii (p > 0), which 
is bounded by the two characteristic heights z equ (q) and 
•2sing(<z)- The boundary z equ (q) agrees with the point of 



9 




equilibrium from Eq. ( 14 ) (see also Fig. [2]) and results 
from root-finding of p(z). The singularity z s ing (g) fol- 
lows from the root of the denominator in Eq. ( |33[ ) and is 
given by 



^sing 



(?) = 



V 



(34) 



a quantity already shown in Fig. [2] The more the ratio q 
of the two masses differs from unity, the more extended 
the interval between z cqu and z s ing becomes. 



Equation (33) follows from 8V e s/dz = with the ef- 
fective potentia For a circular orbit, we also have to 
fulfill dV ca /dp = 0. 



with 



A e (z) = 



1 



2/3 ' 



(37) 



The index "e" designates the case of equal masses. 
Figure 20 shows f e (z) for several M. In the range 
< M < M , with a characteristic mass M , there are 
always two circular orbits at heights z p ± . The numerical 
value of M is given by 



M « 2.598076. 



(38) 



A. Null geodesies 

First, we consider the case in which both masses are 
equal, M := M x = M 2 . The prerequisite dV cB /dp = 
leads to 



[7 + 2,^ = 0. 
dp 



(35) 



Using U(p,z) from Eq. Q and p 2 (z) from Eq. ([33|, the 
expression ( 35 1 yields the root- finding problem 



f c (z) :=[4A c (z)-(l + 2 ) 2 ] 

„ f n M 
-2 2 



M 



M 



\A e {z)Y [A e (z)-z 
M 



[A e (z)Y [A e {z)-z] 



(36) 




M = 0.1 



FIG. 20. (color online) The function f e (z) from Eq. (36 1 for 



several masses M. The heights of the photon orbits are given 
by the roots. Their number depends on the concrete value of 
M. 
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For M > M , there are no intersections with the abscissa 
and no circular orbits. In the limiting case M — M , the 
only root is z p ± — 0, which here does not lie in the range 
of definition. The results of the numerical analysis of 
Eq. ( |36| ) are shown in Fig. 21 
Pp± follow from Eq. 



The corresponding radii 
For each z with < \z\ < 1, 



there is exactly one M G (0, M ) such that z is the height 
of a photon orbit. This is also confirmed by the compar- 
ison of Fig. [2ljb) with Fig. [l^a). 

For unequal masses Mi and M 2 , we find the same equa- 
tion ( 36 ) , but with 



A„ 



Mi 

w 2 



2/3 



1-z 
1+z 



2/3' 



(39) 



instead of A Q (z) from Eq. (37). In this case, we will 



abbreviate the root-finding problem with 
fu(z) = 0. 



(40) 



The solutions will be called z p i with an enumerating in- 



dex i. The corresponding radii are p P i. Figure 22 shows 



/ u (z) for several significant combinations of Mi and M 2 . 
The function depends on both masses, not only on their 
ratio as p(z) from Eq. (33). Because of the symmetry we 



can restrict ourselves to Mi > M 2 at first; thus, we have 
Zcqu < and z s j ng > 0, without loss of generality. 
For f u (z), we find the following limiting values: 



lim f u (z) = -4, lim f u (z) = -oo, 

Z— »Zaing+ 2— S-0 — 

lim f u (z) = oo. 



(41) 



Consequently, there are at least two photon orbits. In 
the interval (— 1;0), the function is monotonic but not 
necessarily for z € (zsing! +!)• Thus, up to four roots can 
exist. Here, the two-dimensional parameter space and 
the complexity of Eq. (40) make exact analysis difficult. 



Significant aspects, however, become clear on the basis 
of the following two examples. 

First, we solve Eq. (40) numerically for variable Mi, 
where M 2 = 1 is fixed (see Fig. 23 1 . There are two photon 



orbits at z p i and z p2 with corresponding radii p p i and 
p p2 . For Mi = 1, the masses are equal and z p i — —z p2 
in agreement with the already discussed case with equal 
masses (see also Fig. 21). 

As a further example, we fix M 2 = 1.8, and Mi re- 
mains variable. The missing monotonicity of f n (z) for 
z € (^sing; +1) for some combinations of Mi and M 2 leads 
to an area of Mi with four photon orbits. Figure[24|shows 
the numerical analysis of Eq. ( 40 1 . 
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FIG. 22. (color online) The function f u (z) from Eq. < |40[ ) for several combinations of the mass values Mi and Mi. The heights 
of the photon orbits are given by the roots. Their number varies between two and four. 



B. Timelike geodesies 



As in the lightlike case, Eq. (33 1 has to be fulfilled. To 
stay on a circular orbit the test particle must have the 
local velocity 



U 



U 



1: 



(42) 



see Appendix [A| Here, we have to use p(z) from Eq. (33) 
and U(p(z), z)hom Eq. ( |4| . For /3 = 1, Eq. |42| con- 
sistently simplifies to Eq. (|35[). As in the lightlike case, 
we distinguish between the cases with equal and unequal 
masses. The function fi(z) is shown in Fig. 25 for some 
values of M := Mi = M^- Since j3 < 1, we can see that 
there are regions of z where no timelike circular orbits are 
possible. The boundaries (f3 = 1) are given by the pho- 
ton orbits, with z p ± calculated in Sec. IV A Because of 
the absence of photon orbits outside the equatorial plane 



for M > M [Eq. (38)], there are also no timelike circular 
orbits outside. 

For the case of unequal masses Mi and M 2 , the velocity 
/3(z) is plotted for several combinations of the mass values 
in Fig. |26| Again, we have an invalid domain between the 
two characteristic values z cqu and z s ing- At the point of 
equilibrium, we have /3(z oqu ) = and p(z equ ) = with 
p(z) according to Eq. (33), i.e., a circular orbit that is 

consistent with Sec. ITTl 



degenerated to a point at 



The appearance of z s - m 
special reason. For z - 



^cqu j 



with f3 — > for z - 
2 s in g +, we have p 



^ s ing-(- has a 
• oo; thus, in 



this limit, there is a circular orbit with infinite radius at 
this height. The velocity decreases with growing radial 
distance, which yields a vanishing velocity in the limiting 
case. These two characteristic heights and the photon 
orbits delimit the domain with possible timclike circular 
orbits. The number of photon orbits can vary between 



two and four (see also Sec. IV A). 
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FIG. 23. (color online) Numerical analysis of Eq. (401 for fixed M2 = 1 and variable Mi. Shown are (a) the heights 2 p i/2 and 
the corresponding radii p p i/2 obtained from Eq. (331 depending on (b) the heights z and (c) the mass Mi. For all Mi, there 
are two photon orbits. 




FIG. 24. (color online) Numerical analysis of Eq. (401 for fixed M2 = 1.8 and variable Mi. Shown are (a) the heights Zpi/2/3/4 
and the corresponding radii Ppi/2/3/4 obtained from Eq. (331 depending on (b) the heights z and (c) the mass Mi. There is a 
range of Mi with four photon orbits (shaded domain). 
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FIG. 25. (color online) Local velocity /3(z) from Eq. (421 for 
a circular orbit at height z. Shown are plots for equal masses 
M := Mi = Mj. The domain with existing photon orbits is 
bounded by the photon orbits discussed in Sec. |IV A| 
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z p+ 
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Zp- 

»M 2 



y 



FIG. 27. (color online) Edge-on view of circular orbits 
for Mi = M2 — 1.5. The photon orbits are located 
at 2 P± « ±0.80324 with p p± w 0.95502. In the plane 
z = 2 s in g = z cqu = 0, there are no photon orbits but an ar- 
bitrary number of timelike orbits (compare with Sees. |III A| 
andlHIBl. 




FIG. 26. (color online) Local velocity /3(z) from Eq. (42 1 for a circular orbit at height z. Shown are plots for several combinations 
of the masses Mi and M%. The domain with existing photon orbits is bounded by the photon orbits discussed in Sec. |IV A] and 
the characteristic heights z cqu and z s in g . 
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y 



FIG. 28. (color online) Edge-on view of circular orbits for 
Mi = 2 and M2 = 1.5. The photon orbits are located at 
z p i » 0.832490 with p pi » 0.80124 and z p2 » 0.35591 with 
p P 2 ~ 1.89345. The plane with z = z S m g ~ 0.14286 marks the 
orbit with infinite radius and vanishing velocity. The point 
on the z axis at z = z cqu « —0.07180 is the circular orbit 
degenerated to a point. 



Mi 



-2p4 



-P 3 Zp2 



~2pl 



"M, 



y 



FIG. 29. (color online) Edge-on view of circular or- 
bits for Mi = 2 and M 2 = 1.8. There are four 
photon orbits with heights z pl w -0.72907, z p2 « 0.08139, 
z P 3 ~ 0.16681, and z P 4 « 0.55673, and radii p p i « 1.00515, 
p p2 » 2.72863, p p3 « 1.81819, and p p4 » 1.31671. The plane 



with 



0.05263 marks the orbit with infinite ra- 



dius and vanishing velocity. The point on the z axis at 



CONCLUSION 



We have analyzed the existence and structure of cir- 
cular orbits in the extreme Reissner-Nordstr0m dihole 
metric for lightlike and timelike test particles and found 
big differences between the different cases. We have also 
shown that such orbits can only exist in heights z between 
the two singularities. 

In the classification of the circular orbits, at first, we 
had to distinguish between two cases. First, we have con- 
sidered a test particle restricted to the equatorial plane of 
the dihole system with equal masses. In this case, the dy- 
namics reduces to an effective two-body problem and the 
structure is similar to a single Reissner-Nordstr0m naked 
singularity investigated by Pugliese et al. jll| . The num- 
ber of photon orbits depends on the concrete mass value: 
There are no orbits for M < M and two for M > M. 
For M = M, we have exactly one photon orbit. Time- 
like circular geodesies are possible for all radii, except 
in the domain between the photon orbits, if such orbits 
exist. Furthermore, we have analyzed the stability be- 
havior in the radial direction. In future work, this could 
be extended to a global stability analysis. 

In the second case, we have focused on circular orbits 
outside the equatorial plane for arbitrary combinations 
of the two mass values. For both lightlike and timelike 
geodesies, their radii are determined by the same expres- 
sion ( 33 ) depending on the height of the orbit and on 
the mass ratio. We have found big differences for equal 
and unequal masses. In the case of equal masses, there 
are two nonequatorial photon orbits for M < M . Their 
heights delimit the range where timelike circular orbits 
are possible. For M > M , there are no photon orbits and 
thus no timelike circular geodesies. For unequal masses, 
in general, there is a forbidden domain between the char- 
acteristic heights z equ and z s i ng . Depending on the com- 



0.02633 is the circular orbit degenerated to a point. bination of the two masses, we find between two and four 



In this case, there are three z domains, where timelike circular 
orbits are possible inside. 



C. Some explicit cases 

We conclude with the presentation of three examples 
that show the three-dimensional structure of lightlike and 
timclikc circular orbits outside the equatorial plane of the 
extreme RN dihole metric. Figures [27ff29] were created 
using the GeodesicViewer by Miiller and Grave 18 . This 
interactive visualization tool numerically integrates the 
geodesic equation with initial conditions with respect to 
a local tetrad system, Eq. The heights z pi of the pho- 



photon orbits. Here, 



and the heights of the 



ton orbits follow from the numerical solutions of Eqs. ( 36 ) 
or ( |40|), respectively, and the corresponding radii p P i from 
Eq. p3[ ). For the timelike orbits, a z value from an al- 
lowed domain is given; the corresponding radius p and 
velocity (3 then follow from Eqs. ( 33 ) and (42 1. 



photon orbits are the boundaries of the ranges with pos- 
sible timelike orbits. In the case of nonequatorial circular 
orbits, further stability analysis would also be of inter- 
est. The extreme Reissner-Nordstr0m dihole metric is 
a static axisymmetric spacetime, and thus the paper of 
Bardeen [T!5] could be helpful in this context. 

In this article, we have studied neutral particles. Be- 
cause of the charges of the two centers of the dihole 
system, trajectories of charged test particles would give 
an even deeper insight into the structure of this space- 
time. Moreover, the Majumdar-Papapetrou spacetimes 
offer the possibility to study geodesies in more complex 
multi black hole configurations. For such discussions, the 
results of this work could be very helpful as reference val- 
ues. Even though such multi black hole spacetimes are 
unlikely to be realized in nature, they could demonstrate 
the diversity of particle orbits in general relativity which 
we will consider in future work. 
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Appendix A: Timelike circular worldline in the 
extreme Reissner-Nordstr0m dihole metric 

The four-velocity u of a timelike particle can be given 
with respect to a local tetrad [see Eq. (10)]. We consider 
a circular worldline (£ = ir/2) with no component in the 
e z direction (x = tt/2). Therefore, the components u p 
and u z have to vanish. The four-velocity then reads 



u = c l ( e (t) + P e ( v )) = l u {p> z ) 9t + 



Pi 



pU{p,z) 



d v (Al) 



with U(p,z) from Eq. Q. The four- velocity u can 
be formally integrated to the circular worldline x^(\) 
parametrized by the affine parameter A. The four- 
velocity u is independent of A because p and z are con- 
stant on the circular worldline. The four- acceleration can 
be calculated according to 



(A2) 



with the Christoffel symbols of the second kind T^g . This 
calculation yields a* = a v = and 



r 
P u 3 



p(l + (3 2 )^L +( 3 2 U 



dp 



1 

TP 



(l- 2 7 2 ) 



dU 

dz 



(A3a) 
(A3b) 



As we are interested in geodesic worldlines, the four- 
acceleration has to vanish. From a z = 0, we obtain 
dU /dz = 0, which leads to the general condition for cir- 
cular orbits, Eq. (15). Equation (A3a) can be solved to 
give 



U 



U- 



- 1 



(A4) 



which is an expression for the local velocity j3 to re- 
main on a circular orbit with radius p. With U{p) from 
Eq. ( 18 ), we also obtain the formula for the local velocity 



in the equatorial plane from Eq. ( 32 ) 



Appendix B: Classical analogon for a massive test 
particle 

In this section, we give a brief overview of the classi- 
cal analogon to the motion of a neutral and massive test 
particle in the extreme RN dihole spacetime. In Fig. [30j 
there is a sketch of this classical system. The masses 
Mi and M 2 have charges Qi an d Qi 01 the same sign. 
Both have the ratio Qi/Mi = 1 (Gaussian and geomet- 
ric units); thus, the gravitational attraction is exactly 
compensated by the electrostatic repulsion. The mass m 
of the test particle is negligible compared to the other 




FIG. 30. (color online) A neutral test particle with mass m is 
attracted by the forces f W and F^ of the two fixed masses 
Mi and M2. For describing a circular orbit, the z components 
of the attracting forces have to compensate each other and the 
p components have to act as centripetal force F z = F~ e p . 



masses. Therefore, Mi and M2 can be considered as 
fixed. 

The test particle is attracted by the forces 



m Mi p 



[p 2 + (z-\yy 

- m M% p 



e a - 



mm{Z - 1 \ e z , (Bla) 

[p 2 + { Z - l) 2 ] 5 

m M 2 (z + 1) 



3 c p 



[ p 2 + {z + iyy [ P 2 + {z + ifY 



-e„ (Bib) 



with the base vectors {ep,e<p,e z } in cylindrical coordi- 
nates. Note that, in this section, bold letters are three- 
dimensional vectors. To describe a circular orbit, the 
z components of the attracting forces must compensate 
each other, which results in the expression 



Mi (z - 1) 
[P 2 + (*-!) 



213 



Mi (z + 1) 
[p 2 + (z + l) 



213 



(B2) 







Equation ( B2 ) is equivalent to the condition dU / dz 
with U{p, z) from Eq. Q. This condition is already 



known from Sec. Ill and can be rearranged to Eq. (15) 



Thus, in this context, the classical system has a similar 
behavior as in the general relativistic case: Circular or- 
bits exist only in the range z £ (—1; +1) and, for z 7^ 0, 
we find an explicit connection between the radii p and the 
heights z of these orbits [see Eq. (33) and Fig. 19 . For 
z = and Mi = M2, the dynamics reduces to an effective 
two-body problem in the equatorial plane. Furthermore, 
to describe a circular orbit, the sum of the p components 
of the attracting forces has to be the centripetal force 



(see Fig. 30). The condition ' +F f 



(2) 



-mv 2 / p, with 



the velocity v on the circular orbit, leads to 



dU 

V 



(B3) 
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FIG. 31. (color online) Classical velocity v(z) from Eq. (B3l for a circular orbit at height z for (a) equal masses M := M\ = M2 
and (b)-(d) several combinations of M\ and M2 with Mi 7^ M-2- As in the relativistic case, the do mai n between the characteristic 
heights 2 G qu and 2 s ing is forbidden. In (b)-(d), the ratios M1/M2 correspond to those from Fig. 26 The speed of light (v = 1) 
is also shown, but represents, in this classical case, no upper boundary of the velocity. 



In this equation, p(z) and U(p(z),z) from Eqs. (33) 
and Q have to be used. Figure [31] shows the velocity 
v(z) on a circular orbit at height z 7M) for several combi- 
nations of Mi and M2 . These functions show a behavior 



qualitatively similar to the relativistic ones of Figs. 25 



and[26j but with the important difference that the speed 
of light, in this classical system, is not an upper boundary 
of the velocity. Thus, in the case of equal masses, there 



are circular orbits for all M and, in the case of unequal 
masses, there are two areas with circular orbits for all 
M 1 and M 2 . 
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